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Abstract 

In this paper we present large-scale correlated calculations of linear optical absorption spectrum 
of oligo-acenes containing up to seven benzene rings. For the calculations we used the Pariser-Parr- 
Pople (P-P-P) Hamiltonian, along with the configuration interaction (CI) technique at various 
levels such as the full CI (FCI), the quadruple CI (QCI) and multi-reference singles-doubles CI 
(MRSDCI). The role of Coulomb parameters used in the P-P-P Hamiltonian was examined by 
considering standard Ohno parameters, as well as a screened set of parameters. A detailed analysis 
of the many-body character of the important excited states contributing to the linear absorption 
has also been performed. The results of our calculations have been compared extensively with the 
theoretical work of other authors, as well as with the experiments. 

PACS numbers: 78.30.Jw, 78.20.Bh, 42.65.-k 
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I. INTRODUCTION 



In last few decades, there has been intensive research in the field of optical properties of 
acenes. -^"^ It is well known that with the increase in the number of benzene rings in the 
acenes, the HOMO to LUMO gap decreases, making the material more conducting. It is 
in fact believed that an infinite linear acene, i.e., the polyacene could be a metal,- and at 
low temperatures, a superconductor.- Thus, higher acenes can prove to be potential candi- 
dates for preparing both optical and electronic devices. Polyacenes generally crystallize in 
well-defined structures, and their crystalline forms have found applications in novel opto- 
electronic devices such as light-emitting field-effect transistors.- Three-dimensional struc- 
tures of polyacenes have also been investigated theoretically for their electronic structure, 
transport and optical properties by several authors.-^ Oligomers of polyacene such as 
naphthalene, anthracene, tetracene, pentacene, etc., are materials which are well-known for 
their interesting optical properties. 10 Because of the high symmetry, they have separate 
optical response to radiation polarized along the conjugation direction, vis-a-vis the radi- 
ation polarized perpendicular to it. Several experimental investigations of linear optical 
properties of polyacenes have been performed over the years. These include linear absorp- 
tion based studies of naphthaleneii^iMiiLIIi^ ant hr acene, ^^^^^^ 
tetracene^amaLataiai pent ^ ne pVMamMMMaaM£L and hexacene.^ Several the- 
oretical studies of the low-lying excited states of these materials have also been performed, 
such as the early LCAO based study by Coulson,— LCAO-MO and perimeter MO approach 
by Pariser,— a free-electron-gas approach by Piatt and Ham et a/.,— ^ CNDO/2 CI approach 
by Hofer et al.f 7 - ab initio multi-reference M0ller-Plesset (MRMP) theory,— density-matrix 
renormalization-group (DMRG) theory using P-P-P model Hamiltonian,^*^ ab initio DFT 
based methodologies,- 4 ^^^ 1 ^ 1 ^ self-consistent field-random phase approximation (SCF- 
RPA) scheme using P-P-P model Hamiltonian,— and CNDO/S2 model parameterization 
technique by Lipari et al? 7 Among the more recent studies using the Pariser-Parr-Pople (P- 
P-P) model Hamiltonian and many-body methodologies, of Ramasesha and coworkers-^ 41 ^ - 
are the foremost. 

Most of the theoretical studies mentioned above concentrate either on a class of excited 
states of polyacenes, or restrict themselves to the study of smaller oligomers. In some earlier 
studies, chiefly because of the lack of computer power at the time, the level of treatment 
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of electron-correlation effects was rather modest by contemporary standards. Additionally, 
to the best of our knowledge, none of the earlier theoretical studies reported calculations 
of optical absorption spectrum of these materials. Therefore, we believe that there is a 
need to perform systematic large-scale many-body calculations of optical properties of these 
materials with the following aims in mind: (i) to compute the linear optical absorption 
spectra of a range of oligoacenes, (ii) to understand their evolution with increasing size of 
the oligomers, (iii) to understand the influence of electron correlation effects on them, and 
finally (iv) to understand the nature of low-lying excited states contributing to the linear 
optics. As it is well-known, in quasi-one-dimensional materials such as conjugated polymers, 
electron-correlation effects have profound influence on their optical properties.— It is with all 
these issues in mind, in the present work, in a systematic manner, we undertake a large-scale 
correlated study of linear optical absorption in oligoacenes of increasing sizes, namely from 
naphthalene to heptacene. The optical absorption spectra have been computed both for the 
radiation polarized along the conjugation direction as well the one polarized perpendicular to 
it. We have used the P-P-P Hamiltonian for the purpose, and utilized various configuration 
interaction (CI) techniques such as the full CI (FCI), the quadruple CI (QCI), and the multi- 
reference singles-doubles CI (MRSDCI) methods. The Cl-based correlation methodology 
used in the present work is sound, and has been used successfully by us in the past to study 
the linear and nonlinear optical properties of various other conjugated polymers.— 
Additionally, we have also examined the issue of the influence of Coulomb parameters on the 
results by performing calculations with two distinct sets of parameters, namely the standard 
Ohno parameters,— and a screened set of parameters meant for phenylene-based conjugated 
polymers,— to describe the P-P-P model Hamiltonian. 

The remainder of this paper is organized as follows. In section [Til we briefly review the 
theoretical methodology adopted in this work. In section DTD we present and discuss our 
results of the optical absorption spectra of various oligoacenes. Finally, in section [TV] we 
summarize our conclusions, and present possible directions for the future research work. 

II. THEORY 

The structures of oligoacenes (C4 n+ 2H 2n +4, n=2, 3, 4, 5, 6, and 7) starting from naph- 
thalene up to heptacene are shown in Figffl Oligomers were assumed to lie in the xy-p\ane 
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with the conjugation direction taken to be along the x-axis. They can be seen as a series of 
benzene rings fused together, along the conjugation direction. An alternative way to look at 
the structure of oligoacenes is to visualize them as two vertically displaced polyene chains 
coupled with each other along the y-axis via hopping, and the Coulomb interactions. From 
this viewpoint, polyacene is a ladder like polymer. The point group symmetry of oligoacenes 
is -D 2 /i, so that the one-photon states belong to the irreducible representations (irreps) B 3u 
or B 2ui while the ground state belongs to the irrep A g . Also, by convention we assign the 
ground state a negative (— ) particle-hole symmetry. Therefore, by dipole selection rules, 
all the optically allowed states must have positive (+) particle-hole symmetry. Thus, the 
states with negative particle-hole symmetry will not be visible in the linear optical spectrum. 
However, we have calculated the energy of IB^ U state, together with the optically allowed 
and 1B^ U states for the sake of comparison with the experimental and other theoret- 
ical results. Clar classified the absorption spectra into three bands, namely, p, a, and j3.— 
In our work, transition to lB^ state from the 1A~ ground state via a short-axis (y-axis) 
polarized photon corresponds to the p band ( 1 L a band of Platt^). Similarly, transition from 
the ground state to 1B% U and 1B£ U via a long-axis (x-axis) polarized photon corresponds to 
the a ( 1 L 6 band of Platt^) and the (5 band of Platt^) bands, respectively. 

The correlated calculations on the oligoacenes were performed using the P-P-P model 
Hamiltonian, which can be written as 



where Hc x and Hc 2 are the one-electron Hamiltonians for the carbon atoms located on the 
upper and the lower polyene chains, respectively. Hc 1 c 2 is the one-electron hopping between 
the two chains, and H ee depicts the electron-electron repulsion. The individual terms can 
now be written as, 



H = H Cl + Hc 2 + H Cl c 2 + H, 



(1) 




(2a) 



(k,k') 




(2b) 



and 
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Figure 1: Schematic drawings of polyacenes considered in this work, namely, (a) naphthalene, (b) 
anthracene, (c) tetracene, (d) pentacene, (e) hexacene, and (f) heptacene 



HciC 2 — t± B k ^. (2c) 
H ee = n i] n H + 2 Vi >i( n i ~ 1 )( n i ~ X ) ( 3 ) 

i " i^j 

In the equation above, k, k' are carbon atoms on the upper polyene chain, /x, v are carbon 
atoms located on the lower polyene chain, while % and j represent all the atoms of the 
oligomer. Symbol (...) implies nearest neighbors, and B it j = J2a( c l.a c j^ + h.c). Matrix 
elements t , and t± depict one-electron hops. As far as the values of the hopping matrix 
elements are concerned, we took t = 2.4 eV for both intracell and intercell hopping, and 
t± = t consistent with the undimerized ground state for polyacene argued by Raghu et al— 
Consequently, the carbon-carbon bond length has been fixed at 1.4 A, and all bond angles 
have been taken to be 120°. At this point it is worthwhile to discuss the issue of the ground 
state geometry of oligoacenes. Experimentally speaking, to the best of our knowledge, the 
available data on the ground state geometry of various polyacenes is for the crystalline 
phase.— Therefore, as far as the ground state geometries of isolated chains are concerned, 

5 



theoretical calculations based upon geometry optimization provide very important input. 
However, the picture which emerges from such calculations is far from clear. Raghu et 
a/.— studied the ground-state geometry of polyacenes using a DMRG based approach and 
concluded that the Peierls' instability in this polymer is conditional, with the gain in the 
electronic energy being proportional to 5 2 (5 is the dimerization amplitude), rather than 
5 which is the case for irans-polyacetylene. In their next paper, Raghu et a/.— concluded 
that the ground state geometry of polyacene consists of a weakly distorted structure with 
undimerized chains, coupled by slightly longer rungs. Therefore, Raghu et al.— used an 
undistorted geometry in their excited state calculations. Houk et al.— , based upon their 
ab initio DFT based calculations, had also predicted a ground state similar to that of 
Raghu et al.— Several workers have suggested that due to Peierls distortion acenes possess 
nonsymmetric geometry,— but on the other hand Cioslowski^ has reported that at the 
correlated level, the symmetric geometry is more stable. Also, Klein and coworkers^ studied 
the distortion of polyacenes within a many-body valence-bond framework, have suggested 
that all the structures (symmetric or nonsymmetric) are close in energy with symmetric one 
possessing the lowest energy. Moreover, Niehaus et al.— have recently studied polyacenes 
using a tight-binding-based Green's-function approach, and reported that for n < 19, the 
symmetric structure is the most stable. Therefore, keeping these uncertainties in mind, 
and the fact that our main aim is to study the influence of electron correlations on the 
optical properties of oligoacenes, we have performed our calculations using the symmetrical 
structure for oligoacenes. This choice also allows us to compare our results directly with 
those of Ramasesha and coworkers,— who also employed a symmetric geometry in their 
excited state of oligoacenes using P-P-P model Hamiltonian. 

The Coulomb interactions are parameterized according to the Ohno relationship,— 

V itj = U/ Kid (l + 0.6117i?y i/2 , (4) 

where, K it j depicts the dielectric constant of the system which can simulate the effects of 
screening, U is the on-site repulsion term, and Rij is the distance in A between the i-th 
carbon and the the j-th carbon. In the present work, we have performed calculations with 
two parameter sets: (a) "standard parameters" with U = 11.13 eV and Kij = 1.0, and (b) 
"screened parameters" with U = 8.0 eV and Kjj = 2.0 (i ^ j) and fi^ = 1, proposed initially 
by Chandross and Mazumdar to study phenyl-based conjugated polymers.— In our earlier 
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studies of phenyl-based conjugated polymers such as the PDPA, PPP, and PPV etc., we 
found that the screened parameters generally provided much better of description of their 
optical properties as compared to the standard ones.— i£2i£M2iS&£i 

The starting point of the correlated calculations for various oligomers were the restricted 
Hartree-Fock (HF) calculations, using the P-P-P Hamiltonian. The many-body effects be- 
yond HF were computed using different levels of the configuration interaction (CI) method, 
namely, full-CI (FCI), quadruples-CI (QCI), and the multi-reference singles-doubles CI 
(MRSDCI). Details of these Cl-based many-body procedures have been presented in our 
earlier works.— From the CI calculations, we obtain the eigenfunctions and eigenval- 
ues corresponding to the correlated ground and excited states of various oligomers. The 
many-body wave functions are used to compute the matrix elements of the dipole operator 
connecting the ground state to various excited states. These dipole matrix elements are in 
turn used to calculate the linear optical absorption spectra of various polyacenes. 

III. CALCULATIONS AND RESULTS 

In this section, first we will briefly discuss the main features of linear optical spectra of 
polyacenes computed within the independent-electron Hiickel model. Next we will present 
and discuss the main results of this work, namely, the correlated linear absorption spectra of 
oligoacenes of increasing sizes, and compare our results with the other available experimental 
and theoretical results. A preliminary description of results for tetracene and pentacene, was 
presented in an earlier work—. However, the results presented here are based upon more 
extensive calculations, and, therefore, they supersede our earlier results-^. 

A. Hiickel Model Results 

Here we briefly discuss the salient features of the linear absorption spectra of oligoacenes 
computed using the tight-binding Hiickel model, and presented in FigfJl 
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Figure 2: Linear optical absorption spectra of (a) naphthalene, (b) anthracene, (c) tetracene, (d) 
pentacene, (e) hexacene, and (f) heptacene. In all the cases a line width of 0.15 eV was assumed. 

For all the oligomers, the first peak (labeled I) corresponds to n — n* excitation described 
by HOMO (H) — > LUMO (L) transition through a ^/-polarized photon leading to the 
excited state of the system. However, with the increasing size of the oligomer, the intensity 
of this H — > L transition decreases, which is understandable because, in the thermodynamic 
limit, with the chosen hopping parameters, the polyacene is metallic. The highest intensity 
peak (labeled II) for each of the oligoacene investigated corresponds to an x-polarized photon 
leading the system to its 1B^ U excited state. For oligomers with n = even, the highest 
intensity peak corresponds to the transitions H —> L + n/2 and H — n/2 — > L, while 
for those with n = odd, the peak II corresponds to transitions H — > L + (n ± l)/2 and 
H — (n ± l)/2 — > L. Moreover, as the lengths of the oligoacenes increase, as expected, the 
spectrum is red shifted due to derealization of particle-hole pair. 
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Table I: The number of reference configurations (N re f) and the total number of configurations 
(Ntotai) involved in the MRSDCI (or FCI or QCI, where indicated) calculations, for different sym- 
metry sub-spaces of various oligoacenes. 



Oligomer 


A 9 

N re f ^total 


B"2u 
N re f N to tal 


B 3u 
N re f Ntotai 


naphthalene 


l a 4936 a 


l a 4794 a 


l a 4816 a 


anthracene 


l a 623576 a 


l a 618478 a 


l a 620928° 


tetracene 


l b 193538 b 


l b 335325 b 


24 c ' d 34788 c ' d 


pentacene 


l b 1002597 6 


l b 1707243 ft 


38 c 130196° 
34 d 126690 d 


hexacene 


66 c 460527 c 
40 d 299141 d 


28 c 191944 c 
32 d 242013 d 


54 c 393248 c 
30 d 252420 d 


heptacene 


45 c 590599 c 
36 d 480032 d 


30 c 415999 c 
22 d 270391 d 


40 c ' d 653476 c ' d 



a FCI method with standard as well as screened parameters, 
b QCI method with standard as well as screened parameters, 
c using standard parameters, 
d using screened parameters. 

B. P-P-P Calculations 

Here we present the results of our correlated calculations of linear absorption on 
oligoacenes using the P-P-P model Hamiltonian. First we present and discuss our results 
for individual acenes, followed by a unified discussion of their spectra. In Table fl] we present 
the number of reference states (N re f) and the dimension of the Hamiltonian matrix (N to tai) 
used in our CI calculations for different symmetry sub-spaces of various oligomers. The fact 
that the calculations presented here are quite large-scale is obvious from Ntotai , which, e.g., 
for pentacene was in excess of one million for the A g and B 2u symmetries. Therefore, we 
are confident that our results take into account the influence of electron-correlation effects 
quite accurately. 

Before discussing the individual oligomers, we would like to summarize the patterns which 
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Table II: Comparison of results of our calculations performed with the standard (Std.) parameters 
and the screened (Scd.) parameters with other experimental and theoretical results for the three 
most important low-lying states. For reference 51 (Kadantsev et ai), 



results quoted with the asterisk (*) correspond to their CISD calculations, while those without it 
are their B3LYP results. 



Excitation energy (eV) 
State Present work Experimental Other theoretical 
Std. Scd. 
para. para. 


Naphthalene (CioHg) 




3.61 


3.22 


3.97-H, 4.0312, 4.0^ 


4.02M, 3.74i£, 3.6012, 4.4451, 4.4951*, 4.2151 4.4652, 4.0912 




4.45 


4.51 


43411, 4.3812, 4.45±£, 4.46H 


4.49il, 4.54#, 4.4612, 4.3551, 5.27-51*, 4.1251,4.8852, 4.6222 


±B+ U 


5.99 


5.30 


5.64H, 5.6212, 5.89I&, 5.95H 


5.94H, 5.84ia, 5.8551, 6.2451*, 5.6951, 5.6222 


Anthracene (C14H10) 




3.25 


2.91 


3.47-n, 3.571a 3<73 aa 


3.72 44 , 3.22 46 , 3.23^«21, 3.8451, 3.9351*, 3.6051, 3.8952 




3.66 


3.55 


3.31-H, 3.3812, 3.42^,25, 3 .4322 


3.65241, 3.60 46 , 3.6811, 3.40^, 3.21-51 4.0551*, 2.9651, 3.6952 




5.34 


4.64 


4.83 11 , 4.8612, 5.2421.25 


5.50H, 5.27^, 5.36H, 4.77a, 5.351, 5. 1451, 5.4951*, 4.9351 


Tetracene (C18H12) 


1^ 


3.22 


3.02 


3.16H, 3.3312, 3.122a 


3.57-H, 3.04l£, 2.9212, 3.4651, 3.5851*, 3.21-51, 3.5252 


1B+ U 


3.16 


2.97 


2.60H^1, 2.71-12, 2.TM, 2.6322 


3.11 44 , 3.051&, 2.8012, 2.4451, 3.2251*, 2.1951, 2.90&* 


^B+ u 


5.01 


4.38 


4.55H, 4.5212, 4.5021, 4.51-22 


5.O9H, 4.86ia, 4.33^, 4.6351, 4.9651*, 4.3851 


Pentacene (C22H14) 




3.17 


2.99 


2.96H, 3.0512, 2.8922, 3.73^ 


3.51-H, 2.991^, 3.2051, 3.3451*, 2.9551, 3.1252 


^B+ u 


2.86 


2.65 


2.1411, 2.2312, 2.1222, 2.28^ 


2.8144, 2.70l£, 1.9051, 2.6451*, 1.6551, 2.37-52 


^B+ u 


4.73 


4.22 


4.01-H, 4.1022, 4.402£ 


4.8011, 4.57-ia, 4.2451, 4.57-51*, 3.9651 


Hexacene (C26H16) 




3.07 


2.77 


2.8C 12 , 2.67i2 


3.3412, 3.01-51, 3.17-51*, 2.7651, 2.87-52 


IB+ U 


2.71 


2.38 


1.90 12 , 1.9li2 


2.1812, 1.5051, 2.2351*, 1.2551, 2.0252 


^B+ u 


4.61 


4.07 


3.9412 


4.05*4, 3.9451, 4.27-51*, 3.6451 


Heptacene (C30H18) 




2.73 


2.35 




2.3551 




2.63 


2.24 


1( 


> 0.9451 


1B+ 


4.48 


3.80 




3.3651 



emerge from the computed optical absorption spectra of acenes ranging from naphthalene 
to heptacene. The absorption spectrum of all the oligoacenes, irrespective of the P-P-P 
parameters used in the calculations, contains following important features: 



1. The first peak in the absorption spectrum for all the oligomers studied corresponds to 
the IB^ U excited state of the system. The most important configuration contributing to 
the many-particle wave function of the state corresponds to the \H — » L) excitation, 
in agreement with Hiickel model calculations. The relative intensity of this feature 
decreases with the increasing size of the oligomer, again in agreement with the Hiickel 
model results. 

2. The second, and the most intense feature in the spectrum corresponds to the 1B% U 
state, obtained by the absorption of a x-polarized photon. The most important con- 
figurations contributing to the wave function of this state, for n = even oligomers are 
excitations \H — > L + n/2) and \H — n/2 — > L), and for n = odd, the excitations 
\H — > L + (n± l)/2) and \H — (n± l)/2 — > L). This aspect of the correlated spectrum 
is also in good agreement with the Hiickel model results. 

3. Another important state, namely 1B^ U state exists for all oligoacenes. Because it 
has the same particle-hole symmetry (— ) as the ground state, in P-P-P (and Hiickel) 
calculations it does not contribute to the absorption spectrum. But many experiments 
report this state as a very weak feature in the absorption spectrum. In our calculations, 
for naphthalene and anthracene, this state occurs at energies lower than the 1B^ U state, 
but for longer oligomers, it is at higher excitation energies than the lB^ state. The 
important configurations contributing to the wave function of this state are the same 
as the ones contributing to 1B% U , except for their opposite relative signs, for all the 
oligomers up to pentacene. But, from hexacene onwards, it is the doubly excited 
configurations which contribute significantly to this state. Thus, it is the electron- 
correlation effects which are responsible for its distinct location in the spectrum as 
compared to the 1B£ U state. 

4. Generally, the spectra computed with the standard parameters are in good qualitative 
agreement with those computed with the screened parameters. Quantitatively speak- 
ing, screened parameter spectra are redshifted as compared to the standard parameter 
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ones, and are in better agreement with the experiments for longer acenes. 

5. Although there are some qualitative similarities between the P-P-P and Hiickel model 
results, quantitatively speaking the optical gaps obtained using the Hiickel model are 
much smaller than their P-P-P and experimental counterparts. 

Having emphasized the general features of our work, next we discuss the results of our 
calculations for individual oligomers, and compare them with the theoretical results of other 
authors, and the experimental ones. Quantitative aspects of our calculations are summarized 
in table HIl which also reports some of the experimental and the theoretical results of other 
authors. Additionally, in Appendix [Aj we present detailed tables for individual oligomers 
containing the excitation energies and many-particle wave functions of important excited 
states, as well as corresponding transition dipoles. 



1. Naphthalene 
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Figure 3: Linear optical absorption spectrum of naphthalene computed with (a) standard parame- 
ters and (b) screened parameters in the P-P-P model Hamiltonian. FCI method was used for the 
purpose. A line width of 0.15 eV was assumed. As a subscript of each feature, its polarization is 
also mentioned. Thus, e.g., Illy implies that peak III is short-axis polarized. 
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In FigJU we present the linear optical absorption spectra of naphthalene computed with 
the FCI method, using both the standard and the screened parameters in the P-P-P Hamil- 
tonian. The wave functions of the excited states contributing to various peaks are presented 
in Tables Hill and [TV] in the appendix. Since these calculations were performed using the 
FCI approach, they are exact within the model chosen, and cannot be improved. Therefore, 
any discrepancy which these results may exhibit with respect to the experiments, is a reflec- 
tion of the limitations of the model or the parameters used, and not that of the correlation 
approach. 

From Fig J3] and Tables HTT1 and II VI it is obvious that the first peak of the spectra calculated 
with the standard and screened parameters are in excellent qualitative and quantitative 
agreement with each other. The second peak computed using the standard parameters is 
due to both the x-polarized and the ^/-polarized components, but it is obvious from the 
Table II I II that x-polarized component dominates. However, with the screened parameters, 
it gets dissociated into two separate features II and III with x-polarized and y-polarized 
components, respectively. Qualitatively, both the spectra are in good agreement with each 
other (see Tables DTD and HVll . Regarding the quantitative aspects, one notices that the 
excitation energies computed with the screened parameters are red shifted as compared to 
those computed with the standard parameters. This is in agreement with results obtained 
for other conjugated polymers as well^j^i^M. 

Comparison of our results for naphthalene, for the invisible (dipole forbidden) 1B^ U state 
and the first two visible features, IB^, and 1B£ U states with the experimental ones, and 
those of other calculations, is presented in Table [III Although, several other experimental 
results on optical absorption in naphthalene exist,— Mdldkd&dlMJJL but here we compare our 
results with the experimental results of Huebner et al— and Aleksandrovsky et a/.,— who 
performed experiments on gas-phase of naphthalene. Since we have performed calculations 
on isolated oligomers, therefore, most appropriate comparison of our results will be with the 
data obtained by gas-phase experiments. Huebner et al.— suggested and IB^ U states 
at 4.45 eV and 5.89 eV respectively, while Aleksandrovsky et al— mentioned these states to 
be at 4.46 eV and 5.95 eV respectively. Thus, our results on the relative ordering of the 1B% U 
and the states in the spectra, using standard parameters are in excellent agreement 
with the results of Huebner et al— and the recent experiment of Aleksandrovsky et al— 
The screened parameters slightly overestimate the experimental results^ 1 ^ for lB^ state 
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and underestimate them for 1B% U state. Excellent agreement between our results and those 
of recent experimental results of Aleksandrovsky et al.— testifies to the essential correctness 
of the P-P-P model to the naphthalene with standard parameters. The dipole forbidden 
lB 3u state has not been discussed by Aleksandrovsky et al. but Huebner et al. suggested 
this feature to be at 4.0 eV which is higher as compared to both of our results. Thus, 
both sets of calculations appear to capture the qualitative features of the spectra quite well, 
however, quantitatively, standard parameter calculations appear to be more accurate upon 
comparison with the experiments.— 1 ^ 

Several theoreticians have studied the low-lying excited states of 
M phthalene<J^fi^l^2ifiifi&SLZMmaa The results obtained from our standard set 
of parameters for the 1B£ U and the 1B% U states are in perfect agreement with the cal- 
culations performed by Pariser et al— and Ramasesha et al.— Also, they agree quite 
well with the recent calculations performed by Rubio and coworkers,— using B3LYP as 
exchange-correlation (xc) functional. The dipole forbidden 1B$ U state is obtained at lower 
energy than calculated by Pariser et alA^ and Rubio and coworkers,— but is in excellent 
agreement with the results of Ramasesha et al.— 

2. Anthracene 

An anthracene molecule contains 14 7r-electrons, which is quite a large system for the 
application of the FCI approach. However, the FCI calculations using the P-P-P model 
have been performed by Ramasesha et alr^ for anthracene, but they did not report its 
linear absorption spectrum. Therefore, in our present work we have computed the linear 
absorption spectrum of the anthracene using the FCI approach, and have also discussed the 
nature of low-lying excited states contributing to the spectrum. 
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Figure 4: Linear optical absorption spectrum of anthracene computed with (a) standard parameters 
and (b) screened parameters. FCI method coupled with the P-P-P model Hamiltonian was used for 
the purpose. A line width of 0.15 eV was assumed. As a subscript of each feature, its polarization 
is also mentioned. 

FigJU presents the linear optical absorption spectra of anthracene computed using the FCI 
method, along with the standard and screened sets of parameters in the P-P-P Hamiltonian. 
The wave functions of the excited states contributing to various peaks are presented in Tables 
IVl and [Vl] in the appendix. 

In our calculations on anthracene using standard parameters, the first peak of the spec- 
trum was identified as 1B% U state at 3.66 eV, while the most intense feature of the spectrum 
corresponding to the 1B^ U state, was obtained at 5.34 eV. The dipole forbidden lB 3u state 
was obtained at 3.25 eV and its wave function contains both singly and doubly excited con- 
figurations. Using screened parameters, we obtained lB 3ul 1B£ U , and 1B£ U states at 2.91eV, 
3.55 eV, and 4.64 eV, respectively. Thus, calculations performed with both sets of P-P-P 
parameters predict the dipole forbidden 1B^ U state at lower energies than the IB^ state. 

Comparison of our theoretical results, to the experimental and theoretical results of other 
authors, is presented in Table [Til Several experimental investigations of the absorption 
spectra of anthracene have been performed over the years. For the l-B^u state, our results 
with both sets of parameters overestimate the experimental results of Klevens and Piatt 
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(3.31 eV)A Biermann et al. (3.38 eV)& Lambert et al. (3.43 eV),^ Dick et al. (3.42 eV)& 
and Man et al. (3.424 eV).— Comparatively speaking, however, the result obtained from 
screened parameters is in slightly better agreement with the experimental results. 

Lambert et alr^ studied the isolated molecules of anthracene using jet spectroscopy, 
while Dick et al.~ observed the low-lying excited states of anthracene in gas-phase. Thus, 
our calculations, which were performed on isolated oligomers, are directly comparable to 
these two experimental results. The most intense feature of the spectra, which corresponds 
to 1B£ U state was obtained at 5.34 eV from standard parameters while screened parameters 
suggested this state to be at 4.64 eV. Our result with standard parameters agrees well with 
the results obtained by Man et a/.,— and Lyons et alr&- They both suggested this state to be 
at 5.24 eV. Our screened parameter result is in good agreement with the result of Klevens 
and Piatt,— and Biermann et al.— for the 1B^ U state, who suggested this state to be at 4.83 
eV and 4.86 eV, respectively. 

Our correlated value of 3.25 eV for the excitation energy of the 1B$ U state obtained using 
the standard parameters is in fairly good agreement with the experimental results of Klevens 
and Piatt (3.47 eV) 11 and Biermann et al. (3.57 eV).— However, it is lower as compared 
to the experimental results of Dick et al. (3.72 eV).— Our screened parameter results (2.91 
eV) obviously significantly underestimate the experimental excitation energy of the 1B% U 
state. Thus, we observe that for some states our standard parameter results agree well with 
the experiments, while for other states screened parameter results are in better agreement 
with the experiments. However, the standard-parameter-based results are in better overall 
agreement with the gas-phase experiments, as compared to those obtained using screened 
parameters. 

Next we discuss and compare our results with the results of calculations performed by 
other theoreticians. Our correlated results obtained by using standard parameters for IB^ U , 
lB^, and 1B£ U states are in perfect agreement with the results obtained by Ramasesha et 
al.^ using FCI methodology with P-P-P model Hamiltonian. Again, this perfect agree- 
ment proves the essential correctness of our calculations. The 1B^ U state obtained using 
standard parameters shows an excellent agreement with the pioneering work of Pariser,— 
who performed the SCI calculations. However, energetically both 1B£ U (dipole allowed) and 
1B$ U (dipole forbidden) states obtained from our calculations are lower than those reported 
by Pariser.— This clearly is due to the superior treatment of electron correlations in our 
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work. Additionally, the results for the l-B^u, and 1B% U states are in excellent agreement 
with results of Hummer et air- and in close agreement with the recent results of Rubio and 
coworkers,— obtained using B3LYP as xc functional. Recently, Kawashima and co-workers^ 
performed calculations on linear absorption in anthracene using ab initio MRMP methodol- 
ogy. The energies they obtained for IB^ state (3.40 eV) and state (4.77 eV) are quite 
close to our screened parameter results, and slightly less than the results obtained from our 
calculations using standard parameters. They mentioned state at 3.23 eV which is 

in perfect agreement with the results which we have obtained using standard parameters, 
while it is higher than the value which we have obtained using screened parameters (2.91 
eV). Few other theoretical results of other authors are also given in the Table HT1 for the sake 
of comparison. 
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Figure 5: Convergence of linear absorption spectrum of tetracene, computed using the MRSDCI 
method (standard parameters) with respect to number of reference configurations (N re f). The 
numbers written in the parenthesis represent N re f included in the corresponding MRSDCI calcu- 
lations. 
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The next member of the polyacene family is tetracene, which contains 18 7r-electrons. 
Because of relatively larger number of 7r-electrons in the system, it is virtually impossible 
to use the FCI method for the present system. Thus, we used QCI technique for computing 
the ground state and the low-lying excited states corresponding to the short-axis polarized 
transitions (f^-type states), while MRSDCI methodology for the excited states which cor- 
respond to the long-axis polarized transitions (-E?3 U -type states). The reason behind using 
different CI techniques is that the QCI method can only be used for single reference states 
such as the ground state (lA g ) and the B 2u type states, while from previous Tables HTT1 HVl 
IVl and fVTl it is obvious that B 3u states are multi-reference states with two dominant con- 
figurations for which QCI method cannot be used. Therefore, for states of B 3u symmetry, 
the MRSDCI method has been used. 
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Figure 6: Linear optical absorption spectrum of tetracene computed with (a) standard parameters, 
and (b) screened parameters. QCI method was used to compute the ground state {IA~ state) and 
B^[ u states, while MRSDCI method was used to compute BJ U states. 

In order to demonstrate the convergence of our MRSDCI calculations, we present the 
spectra in FigJHl calculated with increasing number of reference configurations (N re f). It 
is obvious from the figure that the spectrum computed using twenty four reference config- 
urations is in very good agreement with the spectrum computed using eighteen reference 
configurations, implying that the convergence with respect to N re f has been achieved. We 
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would also like to mention that the accuracy of the MRSDCI method was demonstrated 
in our earlier works where the results obtained using that approach were found to be in 
excellent agreement with the experiments.— 1 ^ 1 ^ 

Tetracene has been studied by several theoreticians^ 1 ^ 1 ^ 1 ^ 1 ^ 1 ^ 1 ^ 1 ^ 1 ^ 1 ^ 1 ^ 1 ^ 1 ^ as well as 
experimentalists.— ii2*2iii2Jli2!i22i22i2i The linear absorption spectra of tetracene, computed us- 
ing standard and screened parameters are presented in FigJBl The excited states together 
with the corresponding energies, transition dipoles, and wave functions are presented in 
Tables ELD and fVTTTl 

The first peak of the correlated spectra was obtained at 3.16 eV and 2.97 eV, using the 
standard parameters and the screened parameters, respectively, and corresponds to the 1B% U 
state of the oligomer. The second peak, which is also the most intense feature in both the 
spectra, corresponds to the 1B^ U state. The dipole forbidden 1B$ U state has been obtained at 
higher energy than the 1B£ U state (c/. Tables [VTT1 and IVlHl) . which is different as compared 
to naphthalene and anthracene, but in agreement with the results of other investigators (cf. 
Table EJ. 

The comparison of our calculated results for 1B^ U , 1B£ U , and 1B£ U excited states with the 
several experimental and other calculated results is present in Table [Til Our standard param- 
eter results generally overestimate the excitation energies when compared to the experiments, 
while the agreement between the screened parameters results and the experiments 11 ^ 1 ^ 1 ^ 1 ^ 
is much better. 

The 1B£ U excitation energies computed using both the standard as well as the screened 
parameters overestimate all the experimental results.— i12i22i22i2L22 Yet our screened param- 
eter value of the excitation energy of 1B% U state (2.97 eV) is in reasonably good agreement 
with the experimental values of Klevens and Piatt (2.60 eV),— Biermann et al. (2.71 eV),— 
Bree and Lyons (2.60 eV)P Birks (2.63 eV)^ and Berlman (2.72 eV)& As far as the 15+ 
state is concerned, our screened parameter value of 4.38 eV is fairly close to the several 
experimental values reported to be near 4.50 eV.— However, our standard parameter 
value of 5.01 eV for the same overestimates the experiments by about 0.5 eV. Our results 
for the dipole forbidden 1B$ U state with both sets of parameters are in good agreement with 
the experiments. As compared to the experiments of Klevens and Piatt 11 - and Berlman^, 
our standard parameter excitation energy for the 1B$ U is slightly higher, while the screened 
parameter value is slightly lower. Thus, we can conclude that in case of tetracene, for 
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all the three excited states discussed, screened parameter based calculations provide better 
agreement with the experiments. 

On comparing our results with other theoretical results we find that IB^ U and 1B^ U states, 
computed using standard parameters are in excellent agreement with the benchmark work of 
Pariser,— while the IB^ U state has been computed at lower value by both of our parameters. 
Ramasesha and co-workers also computed the optical gap of tetracene using P-P-P model 
Hamiltonian coupled with the DMRG technique.— They reported the 1B£ U state at 3.20 
eV, which is in perfect agreement with the 1B% U state computed by us using the standard 
parameters. The results obtained using screened parameters for all the compared states are 
in very good agreement with the results of Ham et al.— and with the ab initio MRMP based 
results of Kawashima et a/.— The results of Rubio and coworkers^ obtained using B3LYP 
as xc functional for the x-polarized states IB^ U and 1B£ U are higher than our results of the 
screened parameters, while is at low energy for the y-polarized 1B£ U state. 

4- Pentacene 

The most thoroughly studied and most widely used member of the polyacene family is 
pentacene. It has been quite famous in recent years due to its use in thin film growth of 
electronic devices.— 122*22 It is among the most promising organic molecular semiconductor 
due to its high charge carrier mobility.-^^ It has been studied by several researchers to 
understand the nature of its low-lying excited states. However, most of the experimentalists 
have mainly reported the energy corresponding to the 1B£ U state.— As per 
our knowledge, experimentally only Klevens and Piatt,— Biermann et a/., 12 and Birks22 
have performed extensive studies of the 1B$ U , lB^, and 1B% U excited states. Additionally, 
Halasinski et al— measured the vibronic transitions of neutral pentacene isolated in Ne, Ar, 
and Kr matrices and reported the data corresponding to the IB^ U , lB^, and 1B£ U excited 
states. 
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Figure 7: Linear optical absorption spectrum of pentacene computed with (a) standard parameters 
and (b) screened parameters. QCI method was used to compute the B^ u states, while MRSDCI 
method was used to compute the B^ u states. 



Theoretically, the low-lying states of pentacene have been studied by several people, 
using various methodologiesT^^&££>££>^ Similar to tetracene, for pentacene also we used 
the QCI method for computing the ground state (lA g ) and the B 2u -tjpe excited states, 
while the B 3u -type excited states were computed using the MRSDCI technique. 

In Figd we present the linear optical absorption spectra computed using both standard 
and screened parameters in the P-P-P model Hamiltonian. The energies and wave functions 
corresponding to the visible features in the spectra as well as state, are presented in 
the Tables [IX] and [X] in the appendix. 

The first peak as usual, corresponds mainly to the lB 2u state, while the second peak is 
a mixture of states corresponding to x-polarized and y-polarized photons. But, as is clear 
from the transition dipoles, the intensity of the peak is mainly due to x-polarized photon 
corresponding to H — * L + 2 + c.c. excitations leading to the l-B^ state (see Tables HXl and 

ED. 

The comparison of our results with other theoretical and experimental results is presented 
in the Table [Til On comparing our results with the experimental results we found that 
on most of the occasions both sets of parameters overestimate the experimental results. 
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However, overall, the results obtained using the screened parameters are in much better 
agreement with the experiments. For example, the screened parameter results on lB^ 
and 1B^ U states are in very good agreement with the experiments of Klevens and Piatt,— 
Biermann et al.,— Birks,— and Halasinski et al.— However, several experimentalists have 
reported that the IB^ U state lies between 1.7 eV to 1.9 eV.— iM.^2.d^M. The difference in our 
calculated results and these experimental results could be due to the bulk effects, as most of 
these experiments were performed on bulk pentacene, thus possibly explaining substantially 
lower values of band gaps. We again emphasize that our calculations were performed on 
single molecule of pentacene, leading to comparatively larger excitation energies. 

Next, we compare our results with the theoretical results of other authors. Our correlated 
results for IB^ and 1B£ U obtained using standard parameters are in very good agreement 
with Pariser— while the energy corresponding to 1-B-^ state is obtained to be lower in our 
standard parameter calculations. Our standard parameter value of lB^ excitation energy 
(2.86 eV) is also in excellent agreement with the value 2.92 eV obtained by Raghu et alM 
from their DMRG calculations. Our results obtained using screened parameters are in very 
good agreement with the results obtained by Ham et a/.— They used free electron molecular 
orbitals (FE MO) theory to compute the energies of the low-lying excited states. From our 
screened parameter calculations, the energy obtained for the IB^ U excited state is in perfect 
agreement with the results of Ham et al.,— while the energies of IB^ and 1B£ U excited states 
are slightly less than that of Ham et a/.— The energies of the x-polarized IB^ U and IB^ U 
states obtained using standard as well as screened parameters agree well with the results for 
these states calculated by Rubio and coworkers,— using xc functional. Recently, Tiago et 
al.— computed the energy gap for solution-phase crystallized (S) structure and vapor-phase 
crystallized (V) structure of pentacene using ah initio pseudo-potential density functional 
method and GW approach. They obtained the energy gap to be 2.2 eV for S structure 
while 1.9 eV for V structure. This low value of energy can be understood due to the bulk 
effects. Thus, their calculation supports our belief that bulk effects affect the H — ► L gap 
significantly. 
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5. Hexacene 



With the increase in the size of oligoacene, it becomes less stable, poorly soluble and more 
reactive.— i22i2L2£i22i2!l Thus, preparation and practical study of higher acenes like hexacene 
is a difficult task— »2i To the best of our knowledge, experimentally hexacene has only 
been studied by Biermann^ and Angliker et al— They measured the absorption spectra of 
hexacene in solution. Angliker et al— have also computed the low-lying excited states by 
performing P-P-P SCF-SCI calculations. Most recently, Grimme et al.,— Heinze et al.,— 
and Houk et a/.— have also calculated the low-lying excited states using the time-dependent 
density functional theory (TDDFT) based approaches. 

In FigJHl we present the linear absorption spectra of hexacene, computed using both 
standard and screened parameters in the P-P-P model Hamiltonian. For inclusion of electron 
correlation effects, the MRSDCI method was employed for all the states involved. The 
excited states, corresponding energies, transition dipoles, and wave functions of the various 
features of the spectra are presented in the Tables ED and IXHI in the appendix. Spectrum 
obtained using standard parameters predicts the 1B^ U state at 2.71 eV, while from screened 
parameters it was obtained at 2.38 eV. In both the cases, the intensity of this peak (peak I) is 
quite small. In the standard parameter spectrum (cf. FigJHh), the second peak corresponding 
to the to the IB^ U state is preceded by a weak shoulder (II y ) which is due to a ^/-polarized 
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photon, and has been identified as state. 
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Figure 8: Linear optical absorption spectrum of hexacene computed with (a) standard parameters 
and (b) screened parameters. MRSDCI method was used to include the electron correlation effects. 

In the screened parameter spectrum, this state is almost degenerate with 1B% U state and 
is part of the main feature II in the spectrum. However, with both sets of parameters, the 
intensity is mainly due to the x-polarized photon (see Tables IXll and IXlB . 

The remarkable feature of the dipole forbidden 1B^ U state for hexacene is that its wave 
function mainly consists of double excitations, with the single excitations making smaller 
contributions. This is quite unlike the smaller acenes in which the 1B$ U state consists 
mainly of singly excited configurations. Thus for this state, contribution of correlation 
effects appears to increase with size of the oligoacene. 

Angliker et ai— have studied the spectra of hexacene experimentally as well as theoreti- 
cally. They dissolved the hexacene in the silicone oil and measured the IB^, lB^, and IB^ U 
excited states at 1.91 eV, 2.67 eV, and 3.94 eV respectively. Theoretically they computed 
these states using P-P-P-SCF-SCI method at 2.18 eV, 3.34 eV, and 4.05 eV respectively. 
Our correlated calculations using standard parameters computed these states at 2.71 eV, 
3.07 eV, and 4.61 eV respectively, while screened parameters found these states to be at 2.38 
eV, 2.77 eV, and 4.07 eV, respectively. Thus, generally our standard parameter calculations 
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overestimate the experimental results of Angliker et a/.— for 1B^ U state. However, with the 
screened parameters, our results are in good agreement with their experimental results^ 
for the IB^ U and 1B^ U excited states. Our screened parameter results for the IB^ and the 
IB^ U states also agree well with their calculated results.— Biermann et al.— observed IB^ U 
at 2.80 eV and IB^ at 1.90 eV. Thus, in very good agreement with our screened parameter 
results for the first DF x-polarized state, while slightly overestimated for the first y-polarized 
state. Again, the trend is clear that the results obtained using the screened parameters are 
in better agreement with the experiments.— 

Recently, Grimme and Parao^ used TDDFT and reported 1B% U and IB^ states at 2.87 
eV and 2.02 eV, which is in close agreement with our screened parameter results. Heinze 
et a/. 5 - calculated the low-lying excited states of polyacenes using TDDFT based coupled 
Kohn-Sham (CKS) methodology. They report IB^, lB^, and 1B% U states to be at 1.25 
eV, 3.64 eV, and 2.76 eV, respectively. Rubio and coworkers— reported these states to be 
at 1.50, 3.94, and 3.01 eV respectively. Except for the 1B£ U state, our results of screened 
parameters calculations are in very good agreement with their results. Additionally, Houk et 
ol.— computed the HOMO to LUMO gap (1-B^ state in our notations) at 1.54 eV, using the 
TDDFT approach. The lower value of 1B£ U excitation energy from the TDDFT calculations 
could be due to the usual problem of underestimating the band gaps, associated with the 
DFT based approaches. 

6. Heptacene 

The final member of the oligoacene family which we have studied is heptacene. It is a big 
molecule with seven diffused benzene rings, and 30 7r-electrons. It is an inherently unstable 
molecule™, whose synthesis has been under controversy— Initially, Clar 94 and Marschalk— 
claimed the preparation of heptacene successfully. Later, Bailey and Liaio reported the 
only other successful synthesis of heptacene.— Just after two years of the report of Bailey 
and Liaio, Clar retracted his earlier claim,— and now it is commonly agreed that heptacene 
cannot be prepared in the pure state.——' Therefore, no experimental data is available till 
date for the low-lying excited states of heptacene. However, the low-lying excited states 
have been studied theoretically by a few researchers using the TDDFT method.——- 
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Figure 9: Linear optical absorption spectrum of heptacene computed with (a) standard parameters 
and (b) screened parameters. MRSDCI method, coupled with the P-P-P model Hamiltonian were 
used for the purpose. 



We computed the linear absorption spectra of heptacene using standard and screened 
parameters in the P-P-P model Hamiltonian, together with the use of MRSDCI technique 
for the inclusion of the correlation effects. The calculated spectra are presented in FigHJ 
while Tables IXIIII and IXIVI in the appendix provide the detailed information regarding the 
excited state properties. We note that, similar to the case of hexacene, the excited state 
wave functions of the and states are fundamentally different in that while the 

1B£ U consists mainly of one-particle excitations, while the IB^ U state is dominated by a 
double excitation, irrespective of the choice of the Coulomb parameters. 

The first feature of FigJHl corresponding to the IB^ state, was obtained at 2.63 eV from 
the standard parameter calculations, and at 2.24 eV on using the screened parameters in 
the P-P-P Hamiltonian. The dipole forbidden state was obtained at 2.73 eV from the 
standard parameters, and at 2.35 eV from the screened parameters. By both the parameters 
IB^ U is at higher energy than as is the case from tetracene onwards. Heinze et 

ai— also calculated the excitation energies for heptacene using the TDDFT based CKS 
methodology.-^ They predicted IB% U state at 0.943 eV, 1B^ U state at 2.349 eV and IB^ U 
state at 3.363 eV. Additionally, Houk et alJ& used TDDFT method and computed the 1B^ U 



26 



state (HOMO to LUMO gap) at 1.24 eV. Our 1B£ U state is higher in comparison to both 
the DFT based results. Lower values of energy for 1B£ U state obtained by DFT based 
methodology could be due to their well known problem of underestimating the band gaps. 
The dipole forbidden 1B$ U states computed by us, using screened parameters is in perfect 
agreement with the results of Heinze et al.,— while the visible 1B% U state with maximum 
intensity is obtained at higher energies by both standard and screened parameters. 

IV. CONCLUSIONS 

In this paper we presented a large-scale correlated study of linear optical absorption, and 
low-lying excited states, of oligoacenes ranging from naphthalene to heptacene, using the 
P-P-P model Hamiltonian. In order to investigate the influence of the Coulomb parameters 
on the computed properties, these calculations were performed with two sets of parameters 
in the P-P-P model: (i) standard Ohno parameters,™^ and (ii) a screened set of parameters.— 
An extensive comparison of our correlated results was made with: (i) the results obtained 
using the single-particle Hiickel model, (ii) the experimental results, and (iii) the theoretical 
results of other authors. Next, in a unified manner, we present the conclusions which can 
be drawn from our calculations. 

Qualitatively speaking, some aspects of the linear optical spectrum computed using the 
Hiickel model, were found to be very similar to the one computed using the P-P-P model, for 
all the oligoacenes investigated. For example, with both sets of calculations, the first peak of 
the spectrum corresponded to the 1B^ U state, while the second, and the most intense, peak of 
the spectrum corresponded to the 1B£ U state for all the acenes. Both sets of calculations also 
predicted diminishing relative intensity of the lB^, with the increasing conjugation length. 
Additionally, irrespective of the parameters used in the P-P-P model, or the size of the acene 
involved, the many-particle wave function of the IB^u state in all the cases consisted mainly 
of the H — ► L excitation, in excellent agreement with the Hiickel model results. However, 
quantitatively speaking, for all the oligoacenes, the optical gaps computed using the Hiickel 
model were much smaller than their P-P-P, and experimental, counterparts. 

As far as the relative position of the dipole forbidden IB^ U excited state is concerned, with 
both the sets of Coulomb parameters the following trends were observed: (i) for oligoacenes 
up to anthracene, 1B^ U occurs below the lB^ state, however, (ii) from tetracene onwards the 
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relative ordering of the two states gets reversed, and state begins to occur below 1B^ U . 
This theoretical result was found to be in good agreement with most of the experiments. 
Additionally, the influence of the electron correlation effects on IB^ U state appears to get 
stronger with the increasing size of the acene involved. This manifests itself in form of 
the trend that the contribution of the singly excited configurations to the many-particle 
wave function of 1B% U decreases with the increasing conjugation length, and is eventually 
overshadowed by the doubly-excited configurations for the larger acenes. 

As mentioned earlier, for all the oligomers studied, irrespective of the Coulomb parameters 
used in the P-P-P Hamiltonian, 1B% U state was found to be the most intense peak of the 
spectrum, in excellent agreement with all the available experiments. Additionally, unlike 
1B$ U state, for all the acenes, the many-particle wave function of this most intense peak 
of the was dominated mainly by the singly-excited configurations. Although, with the 
increasing conjugation length, the relative contribution of the doubly-excited configurations 
to the 1B£ U appears to increase. Therefore, it is conceivable that for conjugation lengths 
longer than the ones considered here, the contributions of the double excitations to IB^ U 
could overshadow those of the single excitations. Yet, for a given conjugation length, the 
electron-correlation effects as judged from the relative contribution of the double excitations, 
appear to be stronger in case of IB^ U state, as compared to the IB^ U state. 

Finally, we discuss the issue pertaining to the choice of Coulomb parameters for the P- 
P-P model, when it comes to describing the linear optical properties of oligoacenes. We 
notice that in an overall sense, with the increasing conjugation length, the results obtained 
using the screened parameters are in a better quantitative agreement with the experiments 
as compared to those obtained using the standard parameters. There were some acenes for 
which both the parameters gave reasonable results, with the screened parameters energies 
slightly lower than the experimental energies, and the standard parameters ones slightly 
higher. Recalling that the screened parameters are generally used to describe the solid- 
state or solvent effects (interchain screening) , the better quantitative agreement obtained 
using the screened parameters for the longer acenes, in our opinion, suggests the increasing 
importance of the solid-state effects on longer acenes. 

In this paper we restricted ourselves to the low-lying excited states of polyacenes which 
contribute to their linear optical properties. However, it will also be interesting to explore 
the nature of their two-photon states, which will contribute to the nonlinear, as well as 
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excited-state absorption, in these materials. Given the inherent anisotropy (short-axis vs. 
long-axis) of these materials, several types of intermediate states will possibly govern their 
nonlinear optical response. At present, studies along these directions are underway in our 
group. 
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Appendix A: EXCITED STATE PROPERTIES 

Here we present the tables summarizing the results of our CI calculations for various 
oligoacenes. The data presented in the tables includes important configurations contributing 
to the many-body wave functions of various excited states, their excitation energies, and 
transition dipoles connecting them to the ground state. The results are presented in separate 
subsections corresponding to each oligoacene, and include calculations performed both with 
the standard parameters, and the screened set of parameters in the P-P-P model. 
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1. Naphthalene 



Table III: Excited states contributing to the linear absorption spectrum of naphthalene computed 
using the FCI method coupled with the standard parameters in the P-P-P model Hamiltonian. 
The table includes many particle wave functions, excitation energies, and dipole matrix elements of 
various states with respect to the ground state. DF corresponds to dipole forbidden state. '+c.c.' 
indicates that the coefficient of charge conjugate of a given configuration has the same sign, while 
'— c.c.' implies that the two coefficients have opposite signs. 



Peak 


State 


E (ev) 


iransition 
Dipole (A) 


tit ttl i * 

Wave functions 


Db 


1B 3u 


3.61 


0.000 


1 TT T i "1 \ i ( f\ /"'OOP' \ 

\H — > L + 1 ) + c.c.(0.6235) 
\H ^ L + 1;H -3^ L) + c.c.(0.1345) 
\H -2 -> L + 3) + c.c.(0.1107) 


T 

1 


1 I ? "1" 


4.45 


0.551 


\H —>■ L) (U.oi 76) 
\H - 1 -» L + 1) (0.3521) 


II 




5.99 


1.638 


|Jf -► L + 1} - c.c.(0.6524) 
|^->L+l;fl"-3-fL) - c.c.(0.1015) 






6.10 


0.739 


|Jf- 1 -» L + 1) (0.8011) 

\H -> L) (0.2958) 
|#-2^L + 2) (0.2952) 


III 




7.78 


1.028 


|#-2^L + 2) (0.8155) 
|Jf- 1 -» L + 1) (0.3002) 


IV 


3£ 3 + 


8.73 


0.281 


|#-2^L + 3) - c.c.(0.4965) 
|ff + - c.c.(0.2280) 
|fT-»L + l;fr-l-»L + 2) - c.c.(0.2260) 
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Table IV: Excited states contributing to the linear absorption spectrum of naphthalene computed 
using the FCI method coupled with the screened parameters in the P-P-P model Hamiltonian. 
The table includes many particle wave functions, excitation energies, and dipole matrix elements of 
various states with respect to the ground state. DF corresponds to dipole forbidden state. '+c.c.' 
indicates that the coefficient of charge conjugate of a given configuration has the same sign, while 
'— c.c.' implies that the two coefficients have opposite signs. 



Peak 


State 


E (eV) 


Transition 
Dipole (A) 


Wave Functions 


DF 


1^ 


3.22 


0.000 


\H -► L + 1} - c.c.(0.5965) 
\H -> L + l-H -?>-> L) + c.c.(0.1564) 


I 




4.51 


0.719 


\H->L) (0.9145) 


II 




5.30 


1.640 


\H -► L + 1} +c.c.(0.6466) 


III 


iBiu 


6.12 


0.842 


\H- 1 -^L + 1) (0.9012) 
\H ^ L + l-H -1 ^ L + 3) - c.c.(0.1334) 


IV 


3£+ 


7.52 


0.844 


\H-2->L + 2) (0.8496) 
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2. Anthracene 



Table V: Excited states contributing to the linear absorption spectrum of anthracene computed 
using the FCI method coupled with the standard parameters in the P-P-P model Hamiltonian. 
The table includes many particle wave functions, excitation energies, and dipole matrix elements of 
various states with respect to the ground state. DF corresponds to dipole forbidden state. '+c.c.' 
indicates that the coefficient of charge conjugate of a given configuration has the same sign, while 
'— c.c.' implies that the two coefficients have opposite signs. 



Peak 


State 


E (eV) 


TV a "n ^ i f i n "n 

J- J. CXJ-lOl 11U11 

Dipole (A) 


Wqvp T^i it) prions 


DF 




3.25 


0.000 


\H — >■ L + 1 } - c.c.(0.5963) 
\H -> L + l-H -?>-> L) - c.c.(0.1329) 
\H -2 -> L + 3) + c.c.(0.1268) 


I 


IBtu 


3.66 


0.728 


\H->L) (0.8894) 
\H- 1 -^L + 1) (0.2097) 


II 




5.34 


2.040 


\H -► L + 1} +c.c.(0.6296) 


III 




6.94 


1.069 


\H -2 -» L + 2) (0.6612) 
|Jf- 1 -> L + 1) (0.4093) 


IV 




7.70 


0.436 


|#-2^L + 3) -c.c.(0.4831) 
\H -> L;H -> L + 2) - c.c.(0.2403) 


V 


7# 2 + 


8.30 


0.769 


|Jf-4 -» L + 4) (0.6955) 
\H-l->L + 5) + c.c.(0.2120) 
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Table VI: Excited states contributing to the linear absorption spectrum of anthracene computed 
using the FCI method coupled with the screened parameters in the P-P-P model Hamiltonian. 
The table includes many particle wave functions, excitation energies, and dipole matrix elements of 
various states with respect to the ground state. DF corresponds to dipole forbidden state. '+c.c.' 
indicates that the coefficient of charge conjugate of a given configuration has the same sign, while 
'— c.c.' implies that the two coefficients have opposite signs. 



Peak 


State 


E (eV) 


Transition 
Dipole (A) 


Wave Functions 


DF 




2.91 


0.000 


| H -» L + 1) - c.c.(0.5670) 
\H ^ L+1;H -3^ L) + c.c.(0.1500) 


I 


. i 


3.55 


0.747 


| H L) (0.8890) 


II 




4.64 


2.019 


| H -» L + 1) + c.c.(0.6215) 


III 




5.93 


0.613 


| H-l^L + 1) (0.7680) 
\H + - c.c.(0.2283) 




3B 2u 


6.01 


0.669 


\H -^L + 4) - c.c. (0.5024) 
\H- 1 -^L + 1) (0.3616) 
\H -2 -» L + 2) (0.2987) 


IV 




6.63 


0.519 


| ii -» L + 5}+c.c. (0.3598) 
| if -► L;H -» L + 2 )+c.c. (0.3170) 
|fT-2->L + 3)-c.c. (0.3161) 


V 


7# 2 + 


7.82 


0.744 


|ii-4 -» L + 4) (0.5739) 
|ii-3^L + 3) (0.3848) 
|ii -» L; if - 1 -> L + 2 )-c.c. (0.2012) 
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3. Tetracene 



Table VII: Excited states contributing to the linear absorption spectrum of tetracene computed 
using the QCI method for A g and B 2u states, and the MRSDCI method for B 3u states, coupled 
with the standard parameters in the P-P-P model Hamiltonian. The table includes many particle 
wave functions, excitation energies, and dipole matrix elements of various states with respect to 
the ground state. DF corresponds to dipole forbidden state. '+c.c.' indicates that the coefficient 
of charge conjugate of a given configuration has the same sign, while '— c.c.' implies that the two 
coefficients have opposite signs. 



Peak 


St cite 


E (eV) 


Transition 
Dipole (A) 


Wave Functions 


DF 


lB 3u 


3.22 


0.000 


|ff -► L + 2} - c.c.(0.5887) 
Iff — 1 — > L 4- 4 ) + c c (0 1395 s ) 
|ff -» L + 2;H - 4 ^ L) - c.c.(0.1289) 


I 


1< 


3.16 


0.806 


\H —> L) (0.8761) 


II 




5.15 


2.392 


\H -^L + 2) + c.c.(0.6198) 


III 




6.35 


0.799 


|ff -2 -» L + 2) (0.5623) 
|ff - 1 -> L + 1) (0.4056) 
|ff-3^L + 3) (0.2211) 
|ff - l^L + 5) + c.c.(0.2066) 


IV 


9# 2 + 


7.90 


0.939 


|ff-3^L + 3) (0.4981) 
|ff-2^L + 6) - c.c.(0.2665) 
|ff-2^L + 2) (0.1925) 
|ff-5^L + 5) (0.1921) 
Iff - 4 -» L + 4) (0.1879) 
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Table VIII: Excited states contributing to the linear absorption spectrum of tetracene computed 
using the QCI method for A g and B 2u states, and the MRSDCI method for B 3u states, coupled 
with the screened parameters in the P-P-P model Hamiltonian. The table includes many particle 
wave functions, excitation energies, and dipole matrix elements of various states with respect to 
the ground state. DF corresponds to dipole forbidden state. '+c.c.' indicates that the coefficient 
of charge conjugate of a given configuration has the same sign, while '— c.c.' implies that the two 
coefficients have opposite signs. 



Peak 


btate 


E (eV) 


Transition 
Dipole (A) 


Wave Functions 


DF 


1^ 


3.02 


0.000 


\H -> L + 2) -c.c.(0.5750) 
\H -2-> L;H -» L + 4) - c.c.(0.1552) 


I 




2.97 


0.799 


\H -» L) (0.8683) 


II 


iB iu 


4.68 


2.356 


\H — »■ L + 2 }+c.c.(0.6269) 


III 


3B 2u 


5.41 


0.671 


| H- 1 -» L+ 1) (0.7657) 
| # -► L + 3)-c.c.(0.2076) 


IV 


,_ I 


6.03 


0.573 


\H—>L]H^L+l }+c.c. (0.4998) 
| #-4^L + l) + c.c.(0.2119) 
|#-1->L+I;fl'->L+1 )+c.c. (0.2004) 






6.10 


0.841 


| H -2->L + 2) (0.7957) 


V 


5£+ 


7.22 


0.413 


| # - 3 -» L + 2 ) -c.c. (0.3768) 
\H -> L;H ^ L+ 1 }-c.c. (0.2426) 
\H ^L;H + 5 }+c.c. (0.2345) 

| if - 1 -» L + 4)+c.c. (0.2177) 






7.37 


0.901 


| F-3^L + 3) (0.7722) 


VI 


165+ 


8.44 


0.650 


| H -b->L + b) (0.6119) 
|^->L+l;fl"-4-^L+l )+c.c. (0.1790) 
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4. Pentacene 



Table IX: Excited states contributing to the linear absorption spectrum of pentacene computed 
using the QCI method for A g and B 2u states, and the MRSDCI method for B 3u states, coupled 
with the standard parameters in the P-P-P model Hamiltonian. The table includes many particle 
wave functions, excitation energies, and dipole matrix elements of various states with respect to 
the ground state. DF corresponds to dipole forbidden state. '+c.c.' indicates that the coefficient 
of charge conjugate of a given configuration has the same sign, while '— c.c.' implies that the two 
coefficients have opposite signs. 



Peak 


State 


E(eV) 


Transition 
DiDole (Al 


Wave Functions 


DF 




3.17 


0.000 


\H -► L + 2}+c.c. (0.5519) 
\H -» L; H -» L + 1 }+c.c. (0.1439) 
\H - 1 — »■ L + 4 )+c c (0 1341) 


I 




2.86 


0.838 


\H ^ L) (0.8568) 


II 




5.01 


2.721 


|if -► L + 2>-c.c. (0.5950) 




3B 2u 


5.16 


0.432 


|iZ"-l-»L + l) (0.5108) 
|if ->■ L + 3)+c.c. (0.3155) 
\H - 2 -► L + 2) (0.2668) 


TTT 

111 




*— i\r\ 

7.09 


0.462 


\H — 1 — ► L + 4 )— c.c. (0.3674) 
\H ^L;H + 5 )+c.c. (0.2329) 
\H ^L;H + 1 )-c.c. (0.1848) 






7.23 


0.413 


|#-1->L+I;fl'->.L + 1 }-c.c. (0.2762) 
\H ^ L;H -1^ L + 3 }-c.c. (0.2266) 
|# -► L + 7}+c.c. (0.2263) 






7.27 


0.668 


|#-3^L + 3) (0.4121) 
|#-6^L + 3)+c.c. (0.2475) 
|#-4^L + 4) (0.2374) 


IV 




7.58 


0.672 


|#-3^L + 3) (0.3294) 
|if - 2 -► L + 2) (0.2682) 
36 \H-5^L + 5) (0.2125) 
|ff-l-»L + 5)-c.c. (0.2060) 



Table X: Excited states contributing to the linear absorption spectrum of pentacene computed 
using the QCI method for A g and B 2u states, and the MRSDCI method for B 3u states, coupled 
with the screened parameters in the P-P-P model Hamiltonian. The table includes many particle 
wave functions, excitation energies, and dipole matrix elements of various states with respect to 
the ground state. DF corresponds to dipole forbidden state. '+c.c.' indicates that the coefficient 
of charge conjugate of a given configuration has the same sign, while '— c.c.' implies that the two 
coefficients have opposite signs. 



Peak 


State 


E (eV) 


Transition 
Dipole (A) 


Wave Functions 


DF 


1^ 


2.99 


0.000 


| H ^ L + 2) + c.c.(0.5117) 
| if ^L;H -> L + 1 ) - c.c.(0.2348) 
if -> L + 2; if - 4 -> L) - c.c. (0.1394) 


I 




2.65 


0.824 


| if -» L) (0.8464) 


II 




4.65 
4.63 


2.642 
0.317 


| if -► L + 2) - c.c.(0.6146) 
| if -► L + 3) + c.c.(0.5123) 
| if -» L) (0.3994) 


III 




6.47 


0.593 


| if - 4 -» L + l)-c.c. (0.4034) 
|if-l^L+l;ff^L + l }+c.c. (0.2226) 
|if -► L; if - 1 -> L + 3 )-c.c. (0.2216) 


IV 




6.77 


0.782 


| if- 3^ L + 3) (0.7460) 
| if -> L;H -> L;H - 3 -» L + 3) (0.2030) 
| ff + - c.c.(0.1413) 
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5. Hexacene 



Table XI: Excited states contributing to the linear absorption spectrum of hexacene computed using 
the MRSDCI method, coupled with the standard parameters in the P-P-P model Hamiltonian. 
The table includes many-particle wave functions, excitation energies, and dipole matrix elements of 
various states with respect to the ground state. DF corresponds to dipole forbidden state. '+c.c.' 
indicates that the coefficient of charge conjugate of a given configuration has the same sign, while 
'— c.c.' implies that the two coefficients have opposite signs. 



Peak 


State 


E (eV) 


Transition 
Dinole ( k) 


Wave Functions 


DF 




3.07 


0.000 


| H -» L-H -» L + 1 ) + c.c.(0.4407) 
| H -► L + 3) +c.c.(0.2729) 
| H ^ L + 7) + c.c.(0.1758) 


I 




2.71 


0.812 


\H -» L) (0.8584) 


TT 
11 


OR+ 
Zn 2u 


A 1 Q 


U. 1 oo 


1 H — ► T, -1- 9 \-i-r r (0 ^04^ 
1 1 r J_j -|- A J^C.C. [U.OVJ^O 1 

| H L;F -1 -> L + 3) + c.c.(0.1810) 
\H ^L;H + - c.c.(0.1720) 


III 


lBt 


4.61 


3.049 


1 # L + 3) - c c f0 5550) 
| H -4-» L + l) +c.c.(0.1653) 


IV 


SB+u 


5.87 


0.764 


| if -► L + 6) +c.c.(0.4457) 
1 77 — 1 — > L 4- 5 \ + r r CO 1 821) 


V 




6.42 


0.549 


>L + 1) + c.c.(0.3513) 
| H-l -^L + 4) +c.c.(0.2390) 
|i7-l^L+l;iJ^L + l} - c.c.(0.2092) 




^B+ u 


6.51 


0.619 


|#-l^L + 5} + c.c.(0.2787) 
\H ^ L + l-H ^ L + 3) + c.c.(0.2784) 


VI 


^B+ u 


6.87 


0.256 


| fl"- 3^ L + 3) (0.4731) 
| H-l ^L + 5) +c.c.(0.2414) 
| H-5^L + 5) (0.2177) 




ZBt u 


6.90 


0.478 


| if -» L; if -» L + 8 ) + c.c.(0.4054) 
|if^L;if-l^L + 6} + c.c.(0.2608) 
| 38-^ L;H -2-> L + 5) + c.c.(0.2090) 


VI 


7B+ U 


7.01 


0.585 


| H-l -^L + 4) +c.c.(0.3445) 



Table XII: Excited states contributing to the linear absorption spectrum of hexacene computed using 
the MRSDCI method, coupled with the screened parameters in the P-P-P model Hamiltonian. The 
table includes many particle wave functions, excitation energies, and dipole matrix elements of 
various states with respect to the ground state. DF corresponds to dipole forbidden state. '+c.c.' 
indicates that the coefficient of charge conjugate of a given configuration has the same sign, while 
'— c.c.' implies that the two coefficients have opposite signs. 



Peak 


State 


E (eV) 


Transition 
Dipole (A) 


Wave Functions 


DF 


1B 3u 


2.77 


0.000 


H -> L;H -> L + 1) + c.c.(0.4943) 
| H -► L + 3) + c.c.(0.1830) 
\H - 1 -» L;H -> L + 2 } - c.c.(0.1819) 


I 


^ u 


2.38 


0.787 


\H -» L) (0.8683) 


II 




3.94 


0.716 


| fl" - 1 -» L + 1 ) (0.7077) 
| H -► L + 2} - c.c.(0.3154) 




i 


4.07 


2.948 


| if -»• L + 3)-c.c.(0.5956) 


III 




5.58 


0.414 


| Lf ->■ L + 6)+c.c.(0.3583) 
| H-l -^L + 5) - c.c.(0.3188) 
| H- 2^ L + 2) (0.2851) 






5.61 


0.971 


| if - 1 ^ L + 4 } + c.c.(0.4188) 
| H -► L + 7)-c.c.(0.2382) 
\H^L;H-l^L + 2 }+c.c.(0.2368) 
|fl'-l->L+l;fl'-fL + 1 )-c.c.(0.2054) 


IV 




5.93 


0.672 


| if- 2^ L + 2) (0.6354) 
| if - 1 -► L + 5)-c.c.(0.2745) 


V 


8B 3u 


6.23 


0.555 


|LL-2^L + 3) + c.c.(0.5539) 
| /L^L + 7)-c.c.(0.1646) 
\H^L;H-l^L + 2 }+c.c.(0.1097) 


VI 


185+ 


7.71 


0.714 


| Lf-3^L + 3) (0.6973) 
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6. Heptacene 



Table XIII: Excited states contributing to the linear absorption spectrum of heptacene computed 
using the MRSDCI method, coupled with the standard parameters in the P-P-P model Hamiltonian. 
The table includes many particle wave functions, excitation energies, and dipole matrix elements of 
various states with respect to the ground state. DF corresponds to dipole forbidden state. '+c.c.' 
indicates that the coefficient of charge conjugate of a given configuration has the same sign, while 
'— c.c.' implies that the two coefficients have opposite signs. 



.r GdK 


otcite 


TP (o\f} 

Hi ^ev J 


ii ciiiDition 
Dipole (A) 


vvdV6 runcooiiD 


Db 


1B 3u 


2.73 


0.000 


1 TT T TT T i 1 \ i / r\ A O C\T*\ 

H — > L; ii — > L + 1) + c.c.(0.4896) 
|#^L + l;iT-2^L} + c.c.(0.1904) 
|ff-»L + l;^-»L + 2) + c.c.(0.1712) 


T 

1 


1B 2u 


C\ /in 

2.63 


0.793 


| .tt — > -L ) (0.8563) 


II 




4.48 


3.251 


| H ^ L + 3) - c.c.(0.5033) 
| —> L; H —>■ L + L)—c.c.{U.2Li\J) 
\H-1 -> L + 5)+c.c.(0.1789) 


III 




5.96 


0.326 


\H ^L;H ^ L + 5 }-c.c.(0.3368) 
\H->L + l-H->L + 3 }+c.c.(0.3312) 
| iT- 1 -» L+ 1) (0.1997) 






6.01 


0.490 


\H->L;H-l->L + 2 }-c.c.(0.4204) 
\H-2^L + 2;H^L+ 1 )-c.c.(0.1781) 


IV 




6.70 


0.624 


| H- 1 ->■ L + 5)+c.c.(0.3258) 
| H -> L;H -1 -> L + 2) - c.c.(0.2972) 
\ H —> L; H —> L + 4) + c.c.(0.2444) 


V 


125+ 


7.09 


0.736 


| F-2^L + 2) (0.3457) 
| F-3^L + 3) (0.3312) 
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Table XIV: Excited states contributing to the linear absorption spectrum of heptacene computed 
using the MRSDCI method, coupled with the screened parameters in the P-P-P model Hamiltonian. 
The table includes many particle wave functions, excitation energies, and dipole matrix elements of 
various states with respect to the ground state. DF corresponds to dipole forbidden state. '+c.c.' 
indicates that the coefficient of charge conjugate of a given configuration has the same sign, while 
'— c.c.' implies that the two coefficients have opposite signs. 



Peak 


State 


E (eV) 


Transition 
Dipole (A) 


Wave Functions 


DF 




2.35 


0.000 


\H -» L;H -» L + l) + c.c.(0.5062) 
| H -1 -» L;H -> L + 2) + c.c.(0.1942) 
\H ^ L + l-H ^ L + 2) + c.c.(0.1648) 


I 


1B 2u 


2.24 


0.795 


H — > L) (0.8675) 


II 




3.80 


1.900 


\H -» L-H -» L+l) - c.c.(0.4795) 
\H-1^L + 1;H^L+1) - c.c.(0.2362) 
\ H -> L + 3 } - c.c.(0.1781) 






3.87 


0.459 


| H ^L + 2) +c.c.(0.5254) 
| H- 1 -» L + 1) (0.4140) 




IBtu 


3.98 


2.553 


| H ^L + 3) - c.c.(0.5710) 


III 




4.92 


0.459 


| H -2 ^ L + 2) +(0.4342) 
|#-l^L + 4} + c.c.(0.4189) 
| # -»■ L + 6 } + c.c.(0.2618) 




±Bt u 


4.98 


0.799 


\H -> L;H -1 -> L + 2) - c.c.(0.4219) 


IV 


SBt u 


5.34 


0.668 


| H ^L + 7) - c.c. (0.4384) 
| H -► L; F -1 -> L + 2) - c.c.(0.2542) 




7B+ U 


5.36 


0.609 


| fl"- 2^ L + 2) (0.5520) 
| # - 1 -» L + 4)+c.c. (0.3475) 


V 




5.57 


0.867 


\H-l->L + 5) + c.c.(0.4741) 
\H ^ L + 7) - c.c.(0.2927) 
\ H ^ L;H -1 ^ L + 2) - c.c.(0.1582) 


VI 


205+ 


8.39 


0.871 


| -» L + 4) (0.7058) 
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